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Abstract. We develop an enhanced technique for characterizing quantum optical processes 
based on probing unknown quantum processes only with coherent states. Our method 
substantially improves the original proposal [M. Lobino et al., Science 322, 563 (2008)], which 
uses a filtered Glauber-Sudarshan decomposition to determine the effect of the process on an 
arbitrary state. We introduce a new relation between the action of a general quantum process 
on coherent state inputs and its action on an arbitrary quantum state. This relation eliminates 
the need to invoke the Glauber-Sudarshan representation for states; hence it dramatically 
simplifies the task of process identification and removes a potential source of error. The new 
relation also enables straightforward extensions of the method to multi-mode and non-trace- 
preserving processes. We illustrate our formalism with several examples, in which we derive 
analytic representations of several fundamental quantum optical processes in the Fock basis. 
In particular, we introduce photon-number cutoff as a reasonable physical resource limitation 
and address resource vs accuracy trade-off in practical applications. We show that the accuracy 
of process estimation scales inversely with the square root of photon-number cutoff. 



PACS numbers: 03.65.Wj, 42.50.-p, 03.67.-a. 
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1. Introduction 

Assembling a complex quantum optical information processor requires precise knowledge of 
the properties of each of its components, i.e., the ability to predict the effect of the components 
on an arbitrary input state. This gives rise to a quantum version of the famous "black box 
problem", which is addressed by means of quantum process tomography (QPT) [Qll2l[3]|. In 
QPT, a set of probe states is sent into the black box (here an unknown completely-positive, 
linear quantum process £ over the set of bounded operators B(Ji) on a Hilbert space W) and 
the output states are measured. From the effect of the process on the probe states it is possible 
to predict its effect on any other state within the same Hilbert space. 

QPT exploits linearity of quantum process over its density operators. If the effect of the 
process £(pi) is known for a set of density operators {pi}, its effect on any linear combination 
p = Y^PiPi equals £(p) = J2Pi£(Pi)- Therefore, if {p^ forms a spanning set within the 
space C(H) of linear operators over a particular Hilbert space H, knowledge of {£(pi)} is 
sufficient to extract complete information about the quantum process. 

However, practical implementations of QPT become demanding especially for systems 
with large Hilbert spaces. For dim('H) = d, dim (£(?{)) = d 2 , which implies that at least d 2 
unknown operators {£(pi}}, each with d 2 unknown parameters, must be estimated. This 
procedure requires preparation of at least {pi}f =1 states, subjecting each to the unknown 
process £, and determining each element of {£(pi)}f =1 through measurement (each with 
d 2 unknown elements), thereby inferring an overall number of d 4 parameters. Furthermore, 
in order to build up sufficient statistics for reliable estimates of the output states, each 
measurement should be performed many times on multiple copies of the inputs. Thus, a 
large number of realizations and measurements are required for complete tomography of £. 

An additional complication, especially for QPT of quantum optical processes, is 
associated with preparation of the probe states. Typical optical QPT implementations deal 
with systems consisting of one or more dual-rail qubits [@l|5l|6]|, which implies that the probe 
states are highly nonclassical, hence difficult to generate. 

These difficulties have been partially alleviated in the recently proposed scheme of 
"coherent- state quantum process tomography" (csQPT) [7J. This scheme is based on the 
observation that the density operator p of a generic quantum state of every electromagnetic 
mode can be expressed in the Glauber-Sudarshan representation ASH!, 



Jc 

where P p (a) is a quasi-probability distribution referred to as the quantum state's "P-function" 
and integrated over the entire complex plane [TTOll . Linearity hence implies that measuring 



i.e., determining the effect of the unknown process on all coherent states enables a prediction 
of its effect upon any generic state p according to 




(1) 



a)(a\ i y Qe{a) = £(\a)(a\), 



(2) 




(3) 
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The implementation of csQPT is advantageous because (i) coherent states are readily 
available from lasers, (ii) coherent states of different amplitudes and phases can be 
produced without changing the layout of the experimental apparatus, and (iii) output-state 
characterization can be performed using optical homodyne tomography ifTTTl . which obviates 
the need for postselection and provides full information about the process in question. 
Moreover, csQPT has been tested experimentally on simple single-mode processes, such as 
the identity, attenuation, and phase shift operations 0. Furthermore, csQPT has been used 
to characterize quantum memory for light based on electromagnetically-induced transparency 

ma. 

An apparent obstacle to csQPT, however, is that the P function for many nonclassical 
optical states exists only in terms of a highly singular generalized function [fT3l[T4l . A remedy 
therefor is provided by Klauder's theorem 031, which states that any trace-class operator p can 
be approximated, to arbitrary accuracy, by a bounded operator pt 6 B(T-L) whose Glauber- 
Sudarshan function P L is in the Schwartz class lfT6l . so integration © can be performed. The 
Klauder approximation can be constructed by low-pass filtering of the P function, i.e., by 
multiplying its Fourier transform with an appropriate regularizing function equal to 1 over 
a square domain of size L x L and rapidly dropping to zero outside this domain. Ref. 
employs this method to implement csQPT. 

Practical implementation of Klauder's procedure is however complicated, because it 
requires finding the characteristic function of the input state and subsequently its regularized 
P function. This function features high-frequency, high- amplitude oscillations that limit 
the precision in calculating the output state ©. Furthermore, Klauder's approximation is 
ambiguous in the choice of the particular filtering function as well as the cutoff parameter 
L0. 

Here we improve csQPT to overcome the above problems. Specifically, we develop 
a new method for csQPT that eliminates the explicit use of the Glauber- Sudarshan 
representation and thus removes the inherent ambiguity associated with employing Klauder's 
approximation for csQPT. In Sec. I2.1[ we obtain an expression for the process tensor in 
the Fock (photon number) basis that can be directly calculated from the experimental data. 
Using this tensor, the process output for an arbitrary input can be calculated by simple 
matrix multiplication rather than requiring integration and high-frequency cut-offs. In this 
way, transformations between the Fock and Glauber-Sudarshan representations, which were 
necessary in Ref. Q, can be sidestepped. Using our new approach, we easily extend csQPT 
from its restrictive single-mode applicability to multi-mode processes and even to non-trace- 
preserving conditional processes. These extensions are particularly relevant for quantum 
information processing circuits, whose basic components are inherently multi-mode and 
conditional ifTTl . 

Process tomography is successful if, for every input state, the estimate for the process 
output closely approximates the actual process output state, and the worst-case error of this 
estimate, given by a distance between the actual and estimated process outputs, is less than 
a given tolerance. For states over infinite-dimensional Hilbert spaces, this concept of error 
is however not meaningful because the finiteness of sampling implies that the process is 
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necessarily under-sampled, hence cannot be determined with bounded error. Instead we could 
consider the process estimation restricted to a finite-dimensional subset of 3(1-1). This version 
of process tomography can always be successful with a sufficiently large amount of sampling. 

Of particular practical interest is the subspace Bifi,) defined by an energy cut-off, i.e., 
estimating the process without accessing any information about its high-energy behavior. This 
restriction is naturally consistent with our choice to work in the Fock basis, because then 
the resulting process tensor is of finite size and with many practical settings (e.g. quantum- 
information processing with photonic qubits). In Sec. 12.21 we provide process error estimates 
for several input state subsets that extend beyond B{W). 

Many interesting processes are phase symmetric; that is, an optical phase shift of the 
input state results in the same phase shift of the output. This property dramatically simplifies 
the experiment because one needs to collect data only for coherent states whose amplitudes lie 
on the real axis rather than the entire complex plane. This prompts us to discuss, in Sec. |3j how 
to obtain the process tensor for phase- symmetric processes, which we test on the experimental 
data from Ref. [fT2| . Next, in Sec. |4l we illustrate our method by analytically deriving the 
superoperators for certain fundamental quantum optical processes using the Fock basis. The 
paper is concluded in Sec. |5] and is supplemented with two appendices. 

2. Coherent state quantum process tomography 

2.1. Formalism: determining the quantum process matrix 

We study general quantum optical processes 8 acting on quantum states of light and begin 
with the simplest case for which only a single electromagnetic field mode is involved. An 
arbitrary quantum state p can be expressed in the Fock representation as 

oo 

p = E Pmn ( n l • ( 4 ) 

m,n=0 

Subjecting this state to an unknown process 8, and imposing linearity, yields 

oo 

£{P)= E Pmn8% n \j)(k\, (5) 
j,k,m,n=0 

where 

8Jl n :=(j\8(\m}(n\)\k) (6) 

is a rank-^ tensor, hereafter referred to as the "process tensor" (superoperator). Thus, by 
expressing input and output states in the Fock basis, a quantum process can be uniquely 
represented and characterized by its rank-^ tensor, which relates the matrix elements of the 
output and input states according to 

[£(p)] Jk = E £ ™p™ ^ 

where N = NU {0}. 
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Below we show how to estimate process tensor elements £( \m) {n\ ) for m, n over a finite 
domain. Because 



,2 a n a m 



(a\(\m)(n\)\a)=e-^^== (8) 

ym\n\ 

is in the Schwartz class, the Glauber-Sudarshan P representation 

\m)(n\ = / d 2 a P mn (a)\a)(a\ (9) 
Jc 

is guaranteed to exist for any operator |m)(n| (m, n £ No) lfT8l . The P function is 

P mn (a) = (-l) m+n ^=d™d n J 2 (a) (10) 

for d™ := d m /da m and a and its complex conjugate a treated as independent variables, and 
5 2 (a) = 5(Re(a))5(lm(a)). By inserting representation © into Eq. ©, and exploiting 
linearity of the process, we obtain the process tensor 

£™ = / d 2 aP mn (a)(j\ 0£ (a)\k). (11) 
Jc 

This expression can be simplified by using Eq. (flOl) and performing integration by parts: 

8 T= I ^^M=Xdl[e\ a \\j\ Qe {a)\k)] 
Jc ym\n\ 



' fd^[e^ 2 (j\ es (a)\k)] 



(12) 

a=0 



\Jm\n\ 

Thus we have eliminated the need to make use of the Glauber-Sudarshan representation 
for quantum states. The process tensor is found by taking partial derivatives (with respect to 
a and a) of the matrix elements of Qs(a), which are estimated from experimental data and 
evaluated at a = 0. 

The mathematical procedure defined by Eq. (fT2~)) is simpler and computationally faster 
(see Sec. [3]) than employing Eq. (fTTT) with a regularized version of PL tmn (a) replacing the 
tempered distribution P mn (a) described in Refs. [13. Equation (fT2~l) has been used 
to determine the fidelity of quantum teleportation of a single-rail optical qubit based on 
measurements performed on coherent states (see supplementary material in Ref. ["1910 . 

Generalization to the multi-mode case is straightforward. In the M-mode case, let us 
introduce the notation |n) := |ra 1 ,n 2 , . . . ,n M ) (with n £ N^) for multi-mode Fock states 
and | a) := \a±, a 2 , • • • , «m) (with a £ C M ) for multi-mode coherent states. Then the matrix 
elements of the output and input states with respect to the Fock basis are related to one another 
by the rank-^ M tensor 

[£( P )} jk = (j\£(p)\k)= Yl £JT/w, (13) 

m,neN* f 

where 

£p k n :=(j\£(\m)(n\)\k). (14) 
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Similarly to the single-mode case, we employ the Glauber-Sudarshan P representation for the 
multi-mode operator \m) (n|, with the overall P function being a product of the P functions 
for the constituent modes: 

Pmn(cc) = J] d™;dZ5 2 (a s ). (15) 

fJl \/m s !n s ! 

Multiple integration by parts yields 

o e (a)\k) 

(16) 

Q S =0 

where 

ee (a)=S(\a)(a\). (17) 

Equations (fl2l) and (fT6l) complete our coherent-state tomography formalism and show that 
coherent states provide a complete set of probe states for characterizing quantum optical 
processes, insofar as the expression for qs{ol) completely determines the process tensor. 

The above formalism is not restricted to trace-preserving quantum processes. Indeed, 
trace preservation was not required in the derivation of our results. Thus, our method is 
applicable to all quantum optical processes that are mathematically described by completely- 
positive maps, but may be trace-preserving, trace-reducing or even trace-increasing. Trace- 
nonpreserving quantum processes are either conditional processes or part of a larger process 
£ = £i + £ 2 , which is trace-preserving as a whole, but whose components £ x and £ 2 
may increase or decrease the trace, respectively. A conditional process is a process that is 
conditioned on a certain probabilistic event; it may be heralded if the event is observed. One 
of the most notable examples of such a process is a probabilistic conditional-NOT gate (CNOT), 
which forms the basis for the Knill-Laflamme-Milburn linear-optical quantum computing 
scheme ifTTl . Other examples are photon-addition and photon- subtraction processes, whose 
superoperators are derived in Sec. HI 

In experimental csQPT, states Qs(a) are determined using homodyne tomography [fTTTl . 
It is important to remember, however, that this procedure reconstructs a density matrix 
normalized to unity trace: Qe(a) = Qs (a)/Tr [ft? (a)]. When measuring non-trace-preserving 
processes, one must recover the trace information contained in Qe(a). This is done by 
measuring the probability p a (£) = Tr [qs (a)} of the process heralding event for all a's for 
which the measurements are performed. The state to be used in Eqs. (fT2l) and (fT6l) in place of 
qs{oc) is then Q £ {a)Ti [gs{a)]. 

An interesting feature of Eqs. (fT2l) and (fT6l) is that complete information about a quantum 
optical process is contained in its action on an infinitesimally small compact set of all probe 
coherent states in the immediate vicinity of the vacuum state. From a mathematical point of 
view, this feature can be understood by realizing that, for any j, k G N , the matrix element 



M 



d 2M c*n , 

(" f^y/m s \n a \ 



» s or, 



till 



M 

n 



e la ^(j\Qs(cx)\k) 
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{j\Qe{ct)\k) is an entire function (see | Appendix A[), i.e., a complex- valued function in the 



variables a, at that is holomorphic over the whole complex plane, and so is its product with the 
exponential e' a ' 2 . Hence, each term e' a ' 2 (j\Qe(ct)\k) is infinitely differentiable over the whole 
complex plane and is identical to its Taylor series expansion in any point of C. Moreover, 
Eq. (fT2l) implies that the process tensor is determined by the corresponding Taylor coefficients 
at a = 0. The same conclusion applies to the multi-mode case, in which we deal with entire 
functions on C M . 

2.2. Energy cutoff and estimation of the error of approximation 

As discussed in Sec. 1, the incompleteness of the information acquired in the experiment is 
accommodated in csQPT by evaluating the process tensor over a restricted finite-dimensional 
subspace "H of the Hilbert space H with a fixed maximum number N of photons. The incurred 
expense is that, through this reduced tomography, only approximate information about the 
process will be inferred: for a given input state p, the predicted output is not £(p), but rather 
£(p), where 

( 18 ) 

Tr[pn] 

is the trace-normalized projection of p onto B(H) and 

£{p) = n£(p)n (19) 



is the predicted output of the reconstructed process for input state p. In Eqs. (1181) and (1191 ), II 
is the projection operator onto H. 

If the input state p is outside B(T-L), the process output estimation error \\£{p) — £ (p ) ||i 
(where ||p|| i = Tr a/ p^p is the trace norm) is generally unbounded. However, it is possible to 
bound the error for certain practically important classes of input states and processes. 

For example, all linear-optical processes involving only linear-optical elements 
(interferometers, attenuators, conditional measurements) do not generate additional photons, 
and thus map B{W) onto itself, so £(p) = £{p). For such processes, the error for a particular 
input p can be estimated according to \\£{p) — < ||£|| — p||i, with the superoperator 

norm defined as ||£|| = sup{||£(.B)||i : B £ B{T-L) , < 1} [20J. If the process is known 

to be trace-nonincreasing, we have ||£|| < 1 [BTI so the error is bounded from above by 

\\£(p)-£(p)\U<\\p-p\\x. (20) 

Note that the above result is not sufficient for evaluating the error for a general process, 
because this error is given by the deviation of £{p) from £(p) rather than from £(p) [Fig. [TJ. 

A further error bound can be obtained for the class of trace-preserving processes that 
do not increase the mean energy, acting on a set of input states whose mean energy does not 
exceed a certain value ||22||. We illustrate this for a single optical mode a with frequency to and 
Hamilton operator H = u(a) a +1/2) whose eigenvalues are denoted by h n = (n + l/2)u. 
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Figure 1. Errors associated with photon number cutoff. Restricting B(H) to B(H) results in 
approximation p of the input state p. If the error of this approximation \\p — p\\i is known, 
the error of the images \\£(p) — £ (p )||i can be estimated according to Eq. ( f2Qb . However, the 
difference between £ (p) and £(p) in the cutoff space remains generally unknown. 



Suppose that the quantum states p of interest satisfy Ti[pH] < U. According to Ref. 11221 . if 
we choose the cutoff dimension dim('H) = N + 1 such that U/h N+ i < 7 for some (small) 
7 > 0, the reconstructed process output errors are bounded from above as 



< 2e, (21) 

1 



where 

6 = 2^7 + 7/(1-7)- (22) 

Conversely, if we want to achieve a certain upper bound e on the error of approximation 
(which corresponds to a lower bound on the desired accuracy of the process characterization), 
we first solve Eq. (|22~1) for 7 = 7(e), and then find the minimum iV 7 G N such that 
U/hjf +i < 7. Any cutoff dimension N + 1 > N 7 is then sufficient for our purpose. For 
7 1, e ~ 2^/7, which yields 

e = 0(1/VN). (23) 



This implies that the error of approximation scales as 1/yN with the cutoff dimension N + 1. 

For example, in order to achieve a 10% error in Eq. (|2~T|) . we need e = 0.05 and thus 
7 ~ 0.0006. For the input mean energy bound corresponding to one photon (U = 3/2u), the 
required cutoff is iV »s {7/7 ~ 250. This calculation shows that the above error estimate is 
extremely conservative. 

3. Phase-invariant processes 



Many practically relevant processes, including the single-mode processes studied in Sec. HI 
exhibit phase invariance. If two input states are identical up to a shift by an optical phase 0, 
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the process outputs for these states differ by the same phase shift: 



-ih<f> 



(24) 



For such processes, it is convenient to express the probe coherent states in polar coordinates: 



| a) 



re 



e |r). Specifically, in these coordinates, we have [|9]| 



\/m\n\ n Am+n 
• r*+i9(n—m) / -j\m+n u S(r) 



(m + n)\ 



and accordingly 



cmn 



y/m\n\ d m+n 



2tt 



d£ 

2^ 



e'" 2 +^— ) 01 &(r,0) | A:) 



{m + n)\ dr m+n 
Hence Eq. ((24)) can be expressed as 

(j\£(\a) (a\)\k)=e^ (j\£(\r) (r\)\k), 

and the superoperator £ [Eq. (|26|) 1 has the following explicit representation: 



r=0 



cmn 



mini 



(m + n)\ dr m+n 



(j\S(\r) (r\)\k) 



8, 



m—j,n—k- 



(25) 



(26) 



(27) 



(28) 



r=0 



In experimental tomography of phase-invariant processes [T2]|, it is sufficient to 
measure the process output for a discrete set of coherent states {|rj)} on the real axis of the 
phase space. The matrix elements of the output states can then be interpolated as polynomial 
functions 



Q 



(j\S(\r) (r\)\k) = Y J Ci{j,k)i 



(29) 



1=0 



where Q is the degree of the polynomial (which depends on the dimension of the truncated 
Hilbert space) and Ci(j, k) are its coefficients. Furthermore, from Eq. (IA.3I ) together with 
Eq. (|28l) . it follows that, for phase-symmetric processes, when j — A; is even or odd, 
(j \ £ (|r) (r|) \k) and its analytic extension to negative values of r are even or odd functions 
of r, respectively. By taking into account the symmetric or antisymmetric property of this 
function, we have additional information to be used in the interpolation procedure; the 
constructed polynomial has to contain only even or odd powers of r, respectively. In this 
way the precision of process estimation from the experimental data is substantially increased. 
With the knowledge of the coefficients Ci(j, k), Eq. (1281) is further simplified to: 



cmn 



y/m\n\ d m+n 
(m + n)\ dr m+n 

I j 7 OO Q 

ymlnl 



si 

s=0 1=0 



5„ 



m—j,n—k 



r=0 



[m + n)l 



~\ C/ ^' k Pm-j,n-k 



s=0 1=0 
l(m+n)}/2 

V m\n\ 

s=0 



si 



Cm+ra-2s(j) k) 



Si 



m—j,n—k 



(30) 
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The last result is significant in that one can obtain the process tensor directly from the 
experimentally reconstructed output states through simple summation. Moreover, if the 
dimension of the truncated Hilbert space is d = N + 1, from Eq. (|30l) it follows that only 
terms of power I < 2N of the interpolation polynomial (|29l ) contribute to the process tensor. 
We have tested this procedure on experimental data [fT2|| and calculated the process tensor in a 
few microseconds, which is a dramatic improvement in comparison to several hours required 
for the original procedure lP7l[T2ll. 

4. Examples: superoperators of important quantum optical processes 

In this section, we illustrate our new method by applying it to some fundamental quantum 
optical processes, whose effects on coherent states are known. Specifically, using Eqs. (fT2)) or 
(fT6l) . we analytically derive corresponding superoperator tensors £^ n in the Fock basis. The 
results are summarized in Tabled! 

4.1. Identity 

For the identity process {£{&), Qs M ( a ) = \ a ) ( a \> me ma trix elements of the output states are 



Inserting these elements into Eq. (PT2l) yields £j£ n = 5 m j5 nk , as expected. 
4.2. Attenuation and lossy channel 

For attenuation of light fields (£ a tt), the process's effect on single-mode coherent states is given 
by 0e M ( a ) = \v a ){v a \> where < rj < 1. The matrix elements in the Fock basis are 



(j\Q£u(<x) \ k ) 



= e 



(31) 




(32) 



From Eq. (fT2~l) . we obtain 



£ 



■mn 





y/m\n\j\k\ 




oo 



(1 -rf) l ai +l a k+l 
Z! 



yjm\n\j\k\ 



QmQuSp 



a,a=0 




(33) 



which depends explicitly on rj. 
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Figure 2. Experimental realizations of (a) photon subtraction and (b) photon addition. The 
process is heralded by single-photon detection events. 



4.3. Photon subtraction and addition 

Photon subtraction is defined as a process that removes a single photon from the light 
field, whereas photon addition adds a single photon. Photon subtraction has been used by 
Ourjoumtsev et al. 11231 to generate optical Schrodinger kittens (coherent superpositions of 
low-amplitude coherent states) from squeezed vacuum states for the purpose of quantum 
information processing. Single-photon- added coherent states can be regarded as the result 
of the most elementary amplification process of classical light fields by a single quantum of 
excitation; being intermediate between single-photon Fock states (fully quantum-mechanical) 
and coherent (classical) ones, these states have been demonstrated to be suited for the study 
of smooth transition between the particle-like and the wavelike behavior of light |[24|. 

Here we discuss idealized single-mode photon subtraction and photon addition. Both 
processes are non-trace-preserving. For example, photon subtraction can be approximately 
realized [23J by a highly-transmissive beam splitter, whose reflected mode is directed to a 
detector and whose transmitted mode constitutes the output, respectively, as illustrated in 
Fig. [2k. Any click in a detector implies extraction of photon(s) from the input mode by 
the beam splitter. As the beam splitter has low reflectivity, here single-photon extraction 
events are more likely than multi-photon events. An approximate experimental realization of 
photon addition is illustrated in Fig. [2b. The input quantum state p enters the signal channel 
of a parametric down-conversion setup. Provided that detector dark counts are neglected, a 
photon detection in the idler mode heralds photon addition to the signal mode, which contains 
the output state of the process. 

The effect of photon subtraction (£ su b) and addition (£ a dd) on coherent states is given by 
Qs^ h {o) = a|a)(a|a''' and Qe. M (ct) = a)\a){a\a, respectively, where a and a) are the photon 
annihilation and photon creation operators of a single mode, respectively. The matrix elements 
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of the output states in the Fock basis are 



(j| fesub (a)|fc) = e-H 2 , (34) 



,12 /— cr a 



01 l*> = e Ha| V^J 7^ = = • 05) 

VU - 1)!(A;- 1)! 



The process tensor is found to be 



en,,, = i a/(j + 1)(A; + l)5 m ,j+i5„,fc+i, for photon subtraction, 
J '* 1 \ffij8m,j-ifin,k-i, for photon addition, 



where we have employed Eq. (fT2)) . 
4.4. Schrddinger cat generation 



The unitary evolution according to £/Keir(x) = exp — i\ (a^a) 2 for \ = 7r /2, if applied to 
coherent states, generates Schrodinger cat states (hereafter denoted as £ cat ) Il25ll26ll 

Q&Jpt) = ^Kerr(|) |«)(«| ^Kerr(|) 

= i(|a)+i|-a))«a|-i(-a|), (37) 

with matrix elements 

01 |fc> = 2/p[ [i + (-iy +k + K-iy - ^(-l) k } . (38) 

The superoperator tensor for this non-Gaussian unitary process obtained via Eq. (fT2l) is 

£$r = e- i *V a -#) 6mj 5 nk . (39) 
Interestingly, this process does not change the total particle number of any input state. 

4.5. Beam splitter 

Now let us consider the beam splitter as an example of a two-mode process. The unitary beam 
splitter transformation is given by [|27ll 

B(Q) = ef^&i-^a^ (40) 

where 6 is the parameter identifying how the beam splitter transmits or reflects beams. 
Specifically, its action on coherent state inputs |ai) and \a 2 ) is given as 

Q £B (a>i,a 2 ) =S B (\ai,a 2 )(ai,a 2 \) 

=&(&)(\a u a 2 }(a u a 2 \)B(Q) 

—\Ta\ — Ra 2 , Ra\ + Ta 2 ) (Ta.\ — Ra 2 , Ra\ + Ta 2 \, (41) 
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with T = cos(6/2) and R = sin(— 6/2) being the transmissivity and reflectivity, 
respectively. By knowing the effect of the process on two-mode coherent states, we can 
calculate the corresponding tensor using Eq. (fT6l) . which yields 



Q 



X I ^ 2 )y 2 P+ 2 9+i2+fc2-mi-ni 

x ^_^ii+fei-p-gj^ji+fei+mi+ni-2p-2g 

X ^mi+m2,ji+j2^ni+«2,fci+fc2' (42) 

as an explicit function of T and R. 
4.6. Parametric down-conversion 

Another two-mode process of interest is parametric down-conversion (PDC). In PDC, a crystal 
with an appreciably large second-order non-linearity is pumped by a laser field. Each of 
the pump photons can spontaneously decay into a pair of identical (degenerate PDC) or 
nonidentical photons (nondegenerate PDC). Here we consider a nondegenerate PDC process 
£ PDC induced by the transformation ll27l 

S 2 (r) = e Kaifc-fiI4). (43) 

The effect of this unitary process on a two-mode coherent state is given by 

£W( a i? «2) = £pDc(|ai, a 2 ) (ax, a 2 \) 

= S 2 {r)\ax,a 2 ) (ax,a 2 \Sl(r). (44) 



In Appendix B we derive the process tensor in the Fock basis. The result can be expressed 
as: 



£ 



m\m2ii\ri2 



'n 1 !m 1 !m 2 !n 2 ! (tanhr) mi+ni Jfl kl 



jlj2klk2 Y jxWhW- (™i - (nx - k x )\ (coshr)^+ fc 2-ji-fei+2 
x 2F1 (-jx,m 2 + 1; mx - jx + 1; tanh 2 r) 
x 2 Fx (—kx, n 2 + l;nx — kx + 1; tanh 2 r) 

x 0rri2— mij'2— ji ^n 2 -ni,k 2 —ki > (45) 



with 



o \ i , (a) n (/3) n z n 
2 F 1 (a,P;r,z) -■= 1 + 2^—rr ) ( 46 ) 

71=1 V ' / n 

the hypergeometric function, (x) n := T(x + n)/T(x) the Pochhammer symbol and T(-) the 
Gamma function [|28l . 
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Table 1. Process tensor £ ^J 1 for some quantum optical processes. 



W JJC1 clLlUll O 


Qs{oi) 


Pi*r\pp>cc tpncnr CTYlTh 


Identity (S^) 


\a)(a\ 


$m j &nk 


Attenuation (£ a tt) 


1 \ / 1 


lm\n\ JJ 3 + fc (l-'7 2 ) m ~ J X 
V j\k\ {m-j)\ °m-j,n-k 


Photon addition (£ a dd) 


a*\oi){a\a 




Photon subtraction (£ su b) 


a\a){a\w 


^ (j + l)(fc + l)5 m<j+1 5 nt k+i 


Cat generation (£ C at) 


\{\a)+i\-a)) 
x({a\ -i{-a\) 




Beam splitter (£b) 


\Toti — Rot2, Rai + Ta.2) 
x (Tax — Ra2 , Ra± + Ta 2 \ 


x( J1 )( J ' 2 )( fel )( k2 ) 

\ p / Vmi — pi \ q / \ni — ql 

w 2^2p+2Q+j 2 +^2 —"Ij —Til 

X ^" 1 +^' 1 + mi + ni -2p— 2g 

^ ^mi+m2 ,ji +J2 i+n2 ,fci+A;2 


Parametric down- 
conversion (£ poc) 




/ mi !m2 !ni In? ! 
V Jilj^lfcilAal 

(mi-ji)! (rii-fci)! (coshr)J2+'=2-Jl-*l+ii 

x 2F1 (— ii, W2 + 1; mi - ji + 1; tanh 2 r) 
x 2-F1 ( — feii n2 + l;ni — fci + 1; tanh 2 r) 



5. Conclusions 

Coherent states are easily generated probe states for tomography of unknown quantum- 
optical processes. Here, we have presented a new, more efficient data processing technique 
for estimating a quantum process from similar experimental procedure of Ref. Q. The 
original formulation was based on regularization and filtering of the Glauber-Sudarshan 
representations for quantum states, which are cumbersome to implement numerically. 
Furthermore, Ref. introduces additional errors associated with regularization of the P 
function. In contrast, our new method to determine the process superoperator [Eq. (PT2l) or 
Eq. (TT6l) 1 is mathematically simpler, computationally faster and unique up to the choice of 
the energy cutoff. Moreover, we presented straightforward generalizations of coherent state 
quantum process tomography to multi-mode and non-trace-preserving conditional processes. 

We have illustrated the new framework through several examples (summarized in 
Tabled)). We have shown that it is straightforward to derive analytically exact and unique 
closed-form expressions for the superoperators for quantum optical processes whose effect 
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on coherent states is known. For phase-invariant unknown processes, the formula to find the 
process tensor reduces to a simple summation of coefficients of a polynomial obtained from 
the experimentally reconstructed output states via interpolation. 

An interesting consequence implied by our formulation [in particular, Eqs. (fT2l) and (fT6l) 1 
is that complete information about a quantum optical process is entirely captured by its effect 
on a compact set of all coherent states |a) in the immediate vicinity of the vacuum state. 
This is due to the entireness property of the image of processes on coherent states. It thus 
appears sufficient to perform tomography experiments only for a range of coherent states 
whose mean photon number is much smaller than that required for the method of Ref. [7] (see 
the suppl. material therein). However, coherent state quantum process tomography relies on 
the ability to approximately determine all the derivatives of a function which is obtained by 
interpolation from measured experimental data. Minimization of errors associated with this 
calculation imposes a lower bound on the phase space region over which the measurements 
need to be performed. For the time being, we have provided an evaluation of the error in the 
process estimation by introducing a truncation of the Fock space. For the class of processes 
respecting a certain energy constraint (which includes all processes that do not amplify the 
energy), we have determined (i) the cutoff dimension that is sufficient in order to achieve 
a certain degree of approximation accuracy, as well as (ii) the upper bound on the error of 
estimation for a given cutoff dimension. 
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Appendix A. Proof that (j\ Qs(ot) \k) is an entire function 

According to Eq. (fT2l) . by knowing the complex- valued function (j \ Qe{a) \k) (of the variable 
a) for any j and k, one can determine the process tensor E™£ 1 . Here we show that this function 
is an entire function so it can be represented as a power series that converges uniformly on 
any compact domain. 

As a completely -positive quantum operation, £ possesses a Kraus decomposition S(p) = 
Yi =1 K iP k\, where L < dim(^) 2 and K { are some Kraus operators on H (whose explicit 
form is not needed for our purpose). Hence we can rewrite the matrix elements of the output 
state as 

L 

(j\g e (a)\k) = J2(j\K i \a) (a\Kj\k) 
i=i 

L 

i=i 

= (a\£*(\k)(j\)\a), (A.l) 
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where 

L 

: B{H) B{H) ,B^^k\Bk h (A.2) 

8=1 

is the dual or adjoint map [|29l . The complex- valued function (a| A |a) (referred to as Husimi 
function if A is a density operator) — where A is any bounded operator on % — is an entire 
function of the two variables a and a lfT3ll30l . Hence the right hand side of Eq. (IA.1I) implies 
that the function (j\ Qs{oi) \k) is an entire function. By representing the coherent states in 
Eq. (IA.1I) in the Fock basis and using Eq. ©, we obtain 

01 6 e(a) \k) = e-H 2 , (A.3) 

n=0 m=0 

which is a power series of the complex variables a and a, hence convergent everywhere 

mum. 

Appendix B. Process tensor for parametric down-conversion 

To obtain the Fock representation of the PDC process, we first find the matrix elements of the 
output states in the Fock basis: 

<Ji , J2 1 £wOi,«2) \h,k 2 ) = (ji, J2I S 2 (r) \a 1: a 2 ) {a x ,a 2 \ S\{r) \k u k 2 ) 

= 1 x J, (B.l) 

where / := (ax, a 2 \ S 2 (r) \k%, k 2 ) and J := (ai, a 2 \ S 2 (r) \ji,j 2 ). Employing the relations 

00 

Sj(r)|0,0) =—. — V(tanhr) z |U>, (B.2) 
cosh r z — ' 
1=0 

S 2 (r)aiS 2 (r) = ai coshr — a^sinhr, (B.3) 

S 2 (r)a 2 S 2 (r) = a 2 coshr — a\ sinhr, (B.4) 

and the binomial expansion, we have: 

J = <«i, « 2 | St(r)^^S 2 (r)St(r) |0, 0) 

Vfci! y/k 2 l 

A //,,\ , 



1 /fc \ 

>(tanh r) 1 (q-, q 2 \ I ) (at coshr) fcl_p (— a 2 sinh 

K2 /A- \ 

y *J (4 cosh r) fc2 - 9 (-ai sinh r) 9 |Z, Z) . (B.5) 



coshr 

k 2 



X 

q=0 
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Using a x \l) = y/l\/(l-x)\ \l - x) and \l) = y/(l + y)\/l\ \l + y) we obtain: 

fci 



- j= | ^(tanhr)' (^J (coshr) fel ~ p (- sinhr) p 



coshr „,.,._. , n 

fc 2 



x ]T ^ (coshr) fc2 - 9 (-sinhr) 9 e- |ai|2/2 ^ |Q2|2/2 

q=0 V ? / 

x (j + *T g)! . (B.6) 

(/ - q)\(l + k 2 -q-p)\ 

From the symmetry between I and J, and by replacing k\ and fc 2 by ji and j 2 , respectively, 
we also find: 

-1 00 h / ■ 

J= 1 V^-— 



— — V (tanhr) 1 ' V (coshr?' 1- ^- sinhrV 

coshrv(riW^ 

x ^ p 2 ") (coshr)^(-sinhrr e -l°il a /2-|a a | 2 /2 



x a? + *— " a£ + *— " 77 (f,T j2 r V)! rr ■ (B.7) 

1 2 (Z / -u)!(Z / + j 2 -w-'y)! 

The Fock representation of the superoperator for the PDC process is then given by 

pm 1 ra 2 nin2 \ ami am f)m 2 an 2 ( Ja-i\ 2 + \a 2 \ 2 j v 7^ 

fcl 



CK1 ,02=0 



-J^T^ , / T ( kl )(coshr) k ^(-smhrY 

V n 2 !m 2 ! coshrVji&Ifcilfo! j^ Q \P J 

x (tanhr) (ni _ fci+p) , <Wn 2 , fe2 -* 2 

x ( 2 j (coshr) fe2_9 (— sinhr tanhr) 9 

<?=0 V ^ / 

x y f Jl Vcoshr)^-»(-sinhr)»(tanhr) mi -^ + " (mi ~ Jl+j2+ ^ )! 
( J (coshr)- ?2_t, (— sinhr tanhr)^ 



u=0 

32 



x ^m 2 -m 1 ,j 2 -ji 

v=0 

In^.m^.k^. (tanhr) ni+mi 



n 2 !m 2 !fc 2 !j 2 ! (coshr) fc2 +^+ 2 

y-y- & r ) kl+n ~ P ~ U (-smhr)P+" 



(n 2 + p)! (m 2 + u)! 

(m - fci (mi - ji + «)! ^-m^Wi > W 
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which can also be expressed in terms of a product of values of the hypergeometric function 
2 Fi, as given by Eq. (1451) . 
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